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The phase diagram of the d — 3 Hubbard model is calculated as a function of temperature 
and electron density ijii), in the full range of densities between and 2 electrons per site, using 
renormalization-group theory. An antiferromagnetic phase occurs at lower temperatures, at and 
near the half- filling density of (n 4 ) = 1. The antiferromagnetic phase is unstable to hole or electron 
doping of at most 15%, yielding to two distinct "r" phases: for large coupling U/t, one such phase 
occurs between 30-35% hole or electron doping, and for small to intermediate coupling U ft another 
such phase occurs between 10-18% doping. Both r phases are distinguished by non-zero hole or 
electron hopping expectation values at all length scales. Under further doping, the r phases yield 
to hole- or electron-rich disordered phases. We have calculated the specific heat over the entire 
phase diagram. The low-temperature specific heat of the weak-coupling r phase shows a BCS- 
type exponential decay, indicating a gap in the excitation spectrum, and a cusp singularity at the 
phase boundary. The strong-coupling r phase, on the other hand, has characteristics of BEC-type 
superconductivity, including a critical exponent a « — 1, and an additional peak in the specific heat 
above the transition temperature indicating pair formation. In the limit of large Coulomb repulsion, 
the phase diagram of the tj model is recovered. 
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I. INTRODUCTION 

The Hubbard model 1] is the simplest realistic (in that 
it retains particulate dynamics) model of electronic con- 
duction systems. This model should constitute a fair 
description for many real solid-state physics systems and 
a starting-point description for those systems with added 
complexities such as quenched randomness, frustration, 
and/or spatial anisotropy. The first query that comes 
to mind, in the study of either experimental or model 
systems, is on the phase diagram, as a function of phys- 
ical parameters such as temperature and density. Nev- 
ertheless, until recently 0, no estimates, let alone (be it 
approximate) solutions, were ventured on the phase dia- 
gram of the Hubbard model at dimensions greater than 
d = 1, at temperatures greater than T = 0, and densities 
away from half-filling. 

The first approach to a phase diagram problem, in the 
past before the advent of renormalization-group theory 
[3j, had been through a mean-field approximation. How- 
ever, such method is not useful for the Hubbard model, 
since, where the characteristic phenomena occur away 
from half-filling, the off-diagonal term in the Hamiltonian 
plays a determining role, as we shall see below. There is 
no ready way to deal with such a dominant quantum 
mechanical effect using mean-field theory. On the other 
hand, renormalization-group theory, which some time 
ago has excelled over mean-field theory in phase diagram 
studies, is effective. Previous renormalization-group cal- 
culations have concentrated on studying the Hubbard 
model in lower dimensions, at zero temperature, or at 
half-filling: The zero-temperature (ground-state) prop- 
erties were successfully obtained in d = 1, 2, 3.(3. 01 In 
d = 1 at half-filling, the thermodynamic properties were 



accurately calculated for finite temperatures. [|j In cases 
where comparison is possible due to the availability of 
exact results in d = 1, the renormalization-group results 
have proven to be very accurate, coming to within about 
1% of the exact results. H In d = 2 at half filling, 
it was found that nophase transition occurs as a func- 
tion of temperature. 013 This result was later extended 
to other fillings in d = 2 [9|| and confirmed by quantum 
Monte Carlo calculations [l(| • In d — 3 at half filling, an 
antiferromagnetic phase transition as a function of tem- 
perature was obtained. One calculation done in d = 3 
at finite temperature and arbitrary chemical potential 
did not obtain the V phase reported below and in 
Ref.@. 

The physics of the Hubbard model in the limit of large 
Coulomb repulsion is believed to be described by the 
t J model [ill \12\ . Application of renormalization-group 
theory to the entire density range of the t J model at 
finite temperatures in d = 3 has yielded 0, 0], be- 
tween 30-40% vacancies from (n,) = 1, a novel (dubbed 
" t" ) phase in which the electron hopping strength in the 
Hamiltonian renormalizcs to infinity under repeated scale 
changes, while the system remains partially filled. The 
calculated topology of the phase diagram, including near 
the r phase a first-order phase transition that is very nar- 
row (less than 2% jump in the electron density) and an 
antiferromagnetic phase that is unstable to at most 10% 
vacancies from (m) = 1, is indeed reminiscent of exper- 
imental phase diagram determinations with lanthanidc 
oxides (lq . 

While the studies above HSI10EH[l![li have 
used position-space renormalization-group approaches, 
there has recently been a revival of interest in Wilson 
perturbative renormalization-group methods applied to 
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correlated fermion problems. These methods have long 
been known to be successful for one-dimensional sys- 
tems 0, 0] and, in the last few years, for the d = 2 
Hubbard model, they have yielded antiferromagnetic in- 
stabilities near half-fillin g an d superconducting instabil- 
ities at smaller densities [il 111 I2I Ell El EsT] . Because 
of the perturbative nature of these treatments, their pre- 
dictions are strictly valid only in the case of weak cou- 
pling. The position-space renormalization-group method 
presented in this paper appears to work over the entire 
range of coupling strengths, as seen below, and yields 
definite phase diagrams and thermodynamic functions. 

In fact, our approach makes an interesting predic- 
tion for the evolution of the Hubbard phase diagram as 
coupling is increased. We find two distinct r phases, 
one occurring at small to intermediate coupling and the 
other, inclusive of the t J model r phase, occurring at 
strong coupling. From an analysis of their specific heat 
behaviors, we find that the two r phases respectively 
have characteristic properties of a weakly-coupled BCS- 
type and a strongly-coupled BEC-type superconducting 
phase. Since high-T c materials share aspects of both lim- 
its, and are thought to lie in some intermediate coupling 
range [24| , our prediction for the Hubbard phase diagram 
may be directly relevant to the physics of high-T c super- 
conductors. 



grouping into a single lattice summation: 



-W =£ -*£ (4c iCT + ct CT c iCT ) (2) 
(ij) { a 

—U (n l i i n il + Ujfriji) + fi (n* + nj)} 

The interaction constants are trivially related by U — 
Uo/2d,fj, = fj,o/2d, and we have hereby exhibited the 
individual-pair Hamiltonian —f3H(i,j). 



III. RENORMALIZATION-GROUP 
TRANSFORMATION 

A. Exact Formulation in d — 1 

For d = 1 (with lattice sites i = 1,2,3,.. .), the Hub- 
bard Hamiltonian in Eq.JSJ takes the form 



II. THE HUBBARD MODEL 

The Hubbard model is defined by the Hamiltonian 

-pH = -t y £ (cl^ + C^cJ) (1) 
(ij)," 

i i 

with f3 = 1/kT, describing electron conduction on a 
d-dimensional hypercubic lattice. Here cJ CT and Ci a re- 
spectively are creation and annihilation operators, obey- 
ing anticommutation rules, for an electron with spin 
a = | or I at the site i of the lattice; = c ia Ci a 
and rii = + are electron number operators. Each 
lattice site can accommodate up to two electrons with 
opposite spins. The index (ij) denotes summation over 
all nearest-neighbor pairs of sites. The three terms of 
this Hamiltonian respectively incorporate kinetic energy 
(parametrized by the electron hopping strength t), on- 
site Coulomb repulsion (with coefficient Uq > 0), and 
chemical potential /io- It is convenient for our purposes 
to rearrange Eq.Q into an equivalent Hamiltonian by 



-PH^Y^i-pH&i+l)}, (3) 

i 

for which an exact renormalization-group transformation 
can be formulated. In terms of matrix elements, this 
exact transformation is fijl 



(uiM 3 u 5 • ■ • |e 13 H |uiu 3 w 5 • • ■ > = 
^ (uiw 2 u i w i u^WQ ■ ■ ■ \e~ tm \ viw 2 v 3 ,w i v 5 w & • • ■ ), 

W2 ,W4,,Wq,.. . 

(4) 

where itj, Wj, and Wi are state variables for lattice site i. 
These variables range over the set {o, f , J,, by which 
we represent the no electron, a single electron with spin 
up, a single electron with spin down, and doubly occu- 
pied states. Here and below, the quantities referring to 
the renormalized (rescaled) system are denoted with a 
prime. The transformation in Eq.Q eliminates half of 
the degrees of freedom in the system, while exactly pre- 
serving the partition function (Z' = Z). However, the 
transformation cannot be readily implemented, due to 
the non-commutativity of the operators in the Hamilto- 
nian. 
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B. Approximation in d — 1 

The renormalization-group transformation formulated 
in Sec.IIIA is implemented approximately, as follows: 



Tr„ 



,-0H =Tr e Ei{-/3H(i,i+l)} 



=Tr„ 



^2Y™{-0H(i-l,i)-l3H(i,t+l)} 



: Y[ Tr J e { ~ /3ff(l ~ 1 ^ ) " ,3ff(4 ' ,;+1)} 



(5) 



-/3'J?'(-i-l,-i+l) 



~ e i:r on {-£'tf'( j -M+i)} = e -/3'ff'. 

In the two approximate steps, marked by ~ in Eq.JSJ), we 
ignore the non-commutation of operators separated be- 
yond three consecutive sites of the unrenormalized sys- 
tem. Since each of these two steps involves the same 
approximation but in opposite directions, some mutual 
compensation can be expected. The success of this 
approximation at predicting finite-temperature behavior 
has been verified in earlier studies of quantum spin sys- 
tems mm. 

The algebraic content of the renormalization-group 
mapping can be extracted from Eq.© as 



-f3'H'(i,k) _ 



-0H(i,j)-0H(J,k) 



(6) 



where i, j, k are three consecutive sites of the unrenormal- 
ized system. The operators —(3'H'(i, k) and —f3H(i,j) — 
f3H(j, k) act on the space of two-site and three-site states 
respectively, so that, in terms of matrix elements, 



[UiV k \e 



-f3'H'(i,k) 



u i v k) 



E 



(ui Wj Vk \e- 0H ^- 0H ^ \Ui Wj v k ) , (7) 



where Ui,Wj,Vk,Ui,v k are single-site state variables. 
Eq.© indicates the contraction of a 64 x 64 matrix on 
the right into a 16 x 16 matrix on the left. This is greatly 
simplified by the use of two- and three-site basis states 
that block-diagonalize respectively the left and right sides 
of Eq.Q). These basis states are the eigenstates of to- 
tal particle number, total spin magnitude, total spin z- 
component, and parity. We denote the set of 16 two-site 
eigenstates by {|0 P )} and the set of 64 three-site eigen- 
states by {IV^}}; and list them in Tables I and II. Eq.J7|) 
is rewritten as 

{M^' H '^ k) \h) = 



v.w 



{■4>q\UiWjV k )(UiVk\4>p) ■ (8) 

In the above equation, with the eigenstates shown in Ta- 



bles I and II, the largest block in {<fi p \e 



-p'H'{i,k) 



I 0p) is 



2x2 and the largest block in (ip q \e' 



-/3H(i,j)-/3H(j,k) 



\lpq) IS 



4x4. (In previous work , some matrix elements in these 
blocks were incorrectly derived). Eq.® yields eleven in- 
dependent elements for the matrix ((ftple - ^ H \cj>p) of 
the renormalized system. These we label 7 P , as shown in 
Appendix A. The values of the 7 P in terms of the matrix 
elements of the unrenormalized system, dictated by the 
right-hand side of Eq. JSJ) , are also given in Appendix A. 
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TABLE I: The two-site basis states used in the derivation 
of the recursion relations, in Eq. (8). In these basis states, 
e -0 h (i,k) j g diagonal, with the exception of a 2 x 2 block 
involving \(f>a) and \<j>a)- The corresponding particle number 
(n), parity (p), total spin (s), and total spin ^-component 
(m 3 ) quantum numbers are also given. The states 1 03) , \4>s), 
|0ii), l^is), 1015} are obtained by spin reversal from \<f>2), 
\4>4,), l^g), |0i2), |0i4), respectively. 



C. Hamiltonian Closed Form under 
the Renormalization-Group Transformation 

Since eleven interaction strengths can be indepen- 
dently fixed by the eleven 7 P , the Hamiltonian —(3'H 1 
which is embodied in Appendix A has a more general 
form than that of the Hubbard Hamiltonian in Eq.©. 
This generalized form of the pair Hamiltonian is 



t h l — (j hj — (j + t x (hi - a rij — (j + n 



E 



l — £7 h j — <T ) 



+t 2 Hi - a nj \c\ a c ja + c] a Ci^ 

- U {n^riii + Tijiriji) + fi (n l + nj) + J S t ■ Sj 
+ V 2 ninj + V 3 (n^n^rij + mn^riji) 
+ ViTi^rnin^nji + G , 



(9) 



where hi a = 1 — m„ is the hole (vacancy) operator 

and Si = J^a a c'ioSo&Ci&i with s aS the vector of Pauli 
spin matrices, is the spin operator at site i. In general, 
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TABLE II: The three-site basis states used in the deriva- 
tion of the recursion relations, in Eq. (8). In these basis 
states, e -t 3H <- i >i)-f 3H U> k ) i s block-diagonal, with the largest 
blocks being 4x4 (see Table IV). The corresponding parti- 
cle number (n), parity (p), total spin (s), and total spin z- 
component (m s ) quantum numbers are also given. The states 

104-5), |07), 1018-19), |022), |027-3o), 1035 - 38), 1041-42), 

1053- 54), 1057), |0eo-6i), 1063) are obtained by spin reversal 

from |02-3), |06), 1014-15), |02o), 1023-26), |031-34), |039-4o), 

|049-5o), 1055), 1058 - 59), 1062), respectively. 



the Hubbard Hamiltonian, after one renormalization- 
group transformation, maps onto this generalized Hamil- 
tonian, which has a form that stays closed under further 
rcnormalization-group transformations. 

The kinetic energy part of the Hamiltonian in Eq.© 
distinguishes the four types of nearest-neighbor hop- 
ping events: i) vacancy hopping (the to term): a 
vacancy (hole) hopping against a background of single- 



electron occupancy (half- filling) ; ii) pair breaking or 
pair making (the t\ term): doubly occupied and com- 
pletely unoccupied nearest-neighbor sites reverting to 
half-filling, or the reverse process; iii) pair hopping 
(the t2 term): a pair hopping against a background 
of half-filling; iv) vacancy - pair interchange (the 
t x term): doubly occupied and completely unoccupied 
nearest- neighbor sites exchanging positions. 

The generalized Hamiltonian of Eq. © reduces to the 
Hubbard Hamiltonian of Eq.Q for t = t\ = t 2 = t and 
t x = J = V2 = V3 = Vi = G = . The renormalization- 
group flows occur in the 10-dimcnsional interaction space 
of the generalized Hamiltonian; the 3-dimensional inter- 
action space of the Hubbard Hamiltonian contains the 
initial conditions of the renormalization- group flows. 

The matrix elements of the renormalized pair Hamil- 
tonian —j3'H'(i, k) are given in Table III in terms of the 
renormalized interaction constants. Table III allows us to 
solve for the renormalized interaction constants in terms 
of the 7 P given in Appendix A: 
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(10) 



where 



1 
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7s - 76 



ln (767s - 7^) 



(78 -7e) +47o 



, cosh 



-1 I 78+76 
2e" 



This completes the determination of our 
renormalization-group transformation, whose flows 
in the ten-dimensional interaction space (to, ii, *2, txj 
U, /Lt, J, V2, V3, V4) are to be analyzed. (G is an additive 
constant not influencing the flows of the 10 other 
interaction constants. However, for expectation value 
calculations, its derivatives must be included in Ea. (|13|) .l 



D. d=l Renormalization-Group Transformation 

The transformation described above is the re- 
moval (decimation) of every other site in a lin- 
ear array. This decimation produces the mapping 
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TABLE III: Block-diagonal matrix of the renormalized two- 
site Hamiltonian —f3'H'(i,k). The Hamiltonian being invari- 
ant under spin-reversal, the spin-flipped matrix elements are 
not shown. 



of a Hamiltonian with interaction constants K = 
(to, ti,t2, t x , U, fx, J, V2,V3, Vi, G) onto another Hamilto- 
nian with interaction constants 



K' = R(K). 



(11) 



The function R is calculated as follows: 

(1) The matrix elements of —j3H(i,j) — f3H(j,k) are 
determined in the three-site basis {ipq} given in Table II. 
In this basis, this matrix is block-diagonal as shown in 
Table IV, with the largest block being 4x4. 

(2) The above block-diagonal matrix is exponentiated, 
yielding the matrix elements (tp q \e~^ H ^ h: '^~^ H ^'' k '\tpq) 
which enter on the right-hand side of Eq.JSJ. This in 
turn yields the eleven j p (as given in Appendix A). 

(3) Using Eqs.JTDJ, the interaction constants of the 
renormalized Hamiltonian —f3'H'(i,k), namely (io> ^11^21 
t' x , U', fj.', J', Vr,', V 3 ', VI, G') are found. 

The initial conditions, for the iterated renormalization- 
group transformations that constitute the 
renormalization-group flow, are the interaction constants 
of the Hubbard Hamiltonian, Ko = (to = t,t\ = t, t^ — 
t, t x =0,U,n,J= 0, V 2 = 0, V 3 = 0, V 4 = 0, G = 0). 



E. d > 1 Renormalization-Group Transformation 

The Migdal-Kadanoff approximation procedure [53, 
l28j (which has been remarkably effective in problems 
as diverse as lower-critical dimensions for different types 
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-J7 + 5m + jJ + 3V 2 + V 3 


^2*2 


>/>52 


V2t 2 


-2(7 + 6m + 4V 2 + 2V 3 





0'55 


V-55 


-U + 5m + 

ij + 3y 2 + v 3 





■056 


■056 


-U + 5m + 
iJ + 3V 2 + V 3 





0>58 


V'59 


■058 


-2U + 6m + 4V 2 + 2V 3 


V2t 2 


'059 


V2t 2 


-3U + 7m + 6V 2 + 5V 3 + Vi 




4>62 






■064 


■062 


-W +7m + 6V 2 + 

5V 3 + Vi 


■064 


-4/7 +8m + 81/ 2 + 

8V, + 2V 4 



TABLE IV: Diagonal matrix blocks of the unrenormalized 
three-site Hamiltonian —/3H(i,j) — /3H(j,k). The Hamilto- 
nian being invariant under spin-reversal, the spin-flipped ma- 
trix elements are not shown. The additive constant contribu- 
tion 2G, occurring at the diagonal terms, is also not shown. 



6 



of phase transitions; first- and second-order phase tran- 
sitions in q-state Potts models; algebraic order in the 
d = 2 XY model; random-field, random-bond, spin-glass 
systems; quenched-disorder-induced criticality; etc.) is 
used to construct the renormalization-group transforma- 
tion for d > 1. In the <i-dimensional hypcrcubic lattice, 
a subset of the nearest-neighbor interactions are ignored, 
so that a hypercubic lattice (still d-dimensional) is left 
behind, in which each lattice point is connected by two 
consecutive nearest-neighbor segments of the original lat- 
tice. The decimation described above can then be applied 
to the site connecting these two segments of the original 
lattice, yielding the renormalized nearest-neighbor cou- 
plings between the lattice points of the new hypercubic 
lattice. To compensate for the nearest-neighbor interac- 
tions that are ignored, the couplings are multiplied by a 
factor of after decimation, b = 2 being the length 
rescaling factor. Thus, the renormalization-group trans- 
formation of Eq. (|ll|l in the previous section generalizes, 
for d > 1, to 



K' = ^-^(K). 



F. Supporting Results 



(12) 



New global phase diagrams obtained by approximate 
renormalization-group transformations are supported by 
the correct rendition of all of the special cases of the 
system solved. The Hamiltonian in Eq. (JSJ, which is the 
system presently solved by approximate recursion rela- 
tions, reduces in various limits to the Ising, quantum 
XY, and quantum Heisenberg spin systems. Our recur- 
sion relations correctly yield the lower critical-dimensions 
di of the Ising (di — 1), quantum XY (di = 2), and 
quantum Heisenberg (di = 2) spin systems. For the 
quantum XY spin system in d — 2, this approximation 
yields the algebraically ordered Kosterlitz-Thouless low- 
temperature phase. [251, |26( For the quantum Heisenberg 
spin system in d = 3, our recursion relations yield low- 
temperature antiferromagnetically (for J < 0) and ferro- 
magnetically (for J > 0) ordered phases, each separated 
by a second-order transition from the high-temperature 
disordered phase. The antiferromagnetic transition tem- 
perature is thus found to be 1.22 times ^| the ferromag- 
netic transition temperature, a purely quantum mechan- 
ical effect, and to be compared with the value of 1.13 
from series expansion |29t l30|. Furthermore, as purely 
off-diagonal quantum effects, the hopping-induced anti- 
ferromagnetism of the d = 3 Hubbard model is recov- 
ered and the scaling of the antiferromagnetic transition 
temperature is obtained with an excellent quantitative 
agreement, as discussed in Sec.V at Ea. lfR))) and shown 
in Fig. 3. In fact, the scaling of the antiferromagnetic 
transition at strong-coupling (Fig. 3), as well as the re- 
sults quoted above, and the disappearance of the transi- 
tion at zero coupling (Fig. 4), indicate the validity of our 
approximation across the entire strong-to-weak coupling 



range. Finally, the Blume-Emery-Grifhths model is con- 
tained in the Hamiltonian of Eq. © and its global phase 
diagram |3l| is obtained from our recursion relations. All 
of these results strongly support the validity of the global 
calculation here. 



IV. RENORMALIZATION-GROUP ANALYSIS: 
GLOBAL PHASE DIAGRAM AND 
OPERATOR EXPECTATION VALUES 

From the recursion equations determined in the pre- 
ceding section, flows are generated for initial values 
of t 1 U, and /i in the Hubbard Hamiltonian. The 
renormalization-group transformation, which constitutes 
each step of the flow, is effected numerically. Particular 
attention has to be given to the multiplication of small 
amplitudes with large exponentials, which can occur in 
the right-hand side of Eq. J5J) when interaction constants 
become large, causing the computational difficulties en- 
countered in previous work 0. 

Each completely stable fixed point, namely sink of the 
renormalization-group flows, corresponds to a thermo- 
dynamic phase, and the global phase diagram is found 
by identifying the basin of attraction for every sink.|3l| 
The expectation values for the operators occurring in the 
Hamiltonian are obtained from the conjugate recursion 
relations, [13 



np = b d n' a T aj 3 , 



(13) 



with summation over the repeated index a implicit. The 
recursion matrix is 



T, 



OK' 



af3 



dK 







(14) 



where K a is an interaction strength, namely a compo- 
nent in the interaction strength vector K defined before 
Eq. (|11|) ; n a is the expectation value of the operator that 
occurs in the Hamiltonian with coefficient K a . Eq. l|13l) is 
iterated along a trajectory until a phase sink limit. The 
left eigenvector of T a p with eigenvalue b d gives the expec- 
tation values at the phase sink, thereby completing the 
calculation of the expectation values of the initial point 
of the trajectory. 

The observed phase sinks in the calculations for the 
d = 3 Hubbard model — the details of which are shown in 
Table V — have a property in common: at the sink limit, 
t\ renormalizes toward zero. In the limit t\ — > 0, analytic 
expressions are derived to first order in t\ for the matrix 
elements l^jj^e -1311 ^^^ 1311 ^^ on the right-hand side 
of Eq.JSJl. This yields, in the neighborhood of each phase 
sink, analytic renormalization-group equations. The ana- 
lytic equations provide a useful check on the accuracy of 
the numerical calculations, and lead to closed-form ex- 
pressions for limiting values of interaction strengths or 
ratios of limiting values of interaction strengths. 
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Flows that start at the boundaries between phases have 
their own fixed points, distinguished from phase sinks by 
having at least one unstable direction. After narrow- 
ing down onto the boundary and from there following 
a flow to the neighborhood of the unstable fixed point, 
a Newton-Raphson procedure is used to exactly locate 
this unstable fixed point. Analysis at these fixed points 
determines the phase transition properties. The expecta- 
tion values calculated, as described above, at the phase 
boundaries allow us to redraw the phase diagram using 
expectation values n a on the axes as well as t, U, and /i. 

The Hamiltonian of Eq. © is covariant under particle- 
hole symmetry (cJ CT — > Ci a ), which in Hamiltonian space 
takes the form of a mapping K = S(K). The function S 
is given by 

to = — t 2 , t\ = — 1\, t 2 = —to, t x = tx> J = J, 
U = U - 2V 3 - Vi, p = -li + U - 2V 2 - W 3 - Va, 
V 2 = V 2 + 2V 3 + V 4 , V 3 = -V 3 - Vi, V 4 = V 4 . (15) 

The subspace that is invariant under S corresponds to 
systems that are invariant under particle-hole exchange, 
and therefore are at half-filling: (rii) = 1 = (hi). From 
Ea.lfT5)l. this subspace occurs at to = —t 2 , t\ =0, 2\i = 
U - 2V 2 - V 3 , 2V 3 = -V 4 . For the original Hubbard 
Hamiltonian, all points with fio/Uo = 1/2 are mapped 
onto this subspace after the first renormalization-group 
step. 

The Hubbard phase diagrams are plotted in the 
next section, for fixed Uo/t, in terms of 1/t (a tem- 
perature variable) versus [lo/Uo or (rii). Since our 
renormalization-group transformation is also covariant 
under particle-hole symmetry, the phase diagrams are 
duly symmetric about Ho/Uq — 1/2 or (m) = 1. 

V. GLOBAL PHASE DIAGRAM FOR d = 3 

For d = 3 and a range of couplings Uo/t = 5 to 20, 
Figs. 1 show Hubbard phase diagrams in terms of tem- 
perature (1/t) versus chemical potential (ho/Uq). The 
corresponding phase diagrams in temperature (1/t) ver- 
sus electron density (rii) are in Fig. 2. The values of the 
interaction constants for each observed phase sink are 
listed in Table V. The expectation values for each phase 
sink, also listed in Table V, allow us to identify the phases 
as follows: 

Hole-rich disordered (hD) phase: The electron 
density (m) is zero at the sink and, concomitantly, the 
electron densities (n») calculated inside this phase are 
low. 

Near-half-filled disordered (nHD) phase: The 

basin of attraction of nHD occurs at Ho/Uo ^ 1/2. The 
electron density (rii) is 1 at the sink and, concomitantly, 
the electron densities (rii) calculated inside this phase 
are closer to half-filling. Half-filled disordered (HD) 
phase: The sink is for the disordered phase at perfect 
half-filling, (Iq/Uq = 1/2 and (m) = 1. 



Electron-rich disordered (eD) phase: The elec- 
tron density (rii) is 2 at the sink and, concomitantly, the 
electron densities (rii) calculated inside this phase are 
high. 

Antiferromagnetic (AF) phase: The electron den- 
sity (rii) is 1 at the sink and, concomitantly, the elec- 
tron densities (rii) calculated inside this phase are closer 
to half-filling. The expectation value for the nearest- 
neighbor spin-spin correlation is (Si ■ Sj) = \ at the sink. 
Note that the latter two spins are, on the original cu- 
bic lattice, distant spins on the same sublattice; from 
this, antifcrromagnetism, (Si ■ Sj) < when the spins 
are on different sublattices of the original cubic lattice, 
is calculationally obtained throughout this phase. Since 
there is no explicit antiferromagnetic coupling in the ini- 
tial Hubbard Hamiltonian, the antiferromagnetic phase 
is completely a quantum mechanical effect resulting from 
the kinetic energy term. In fact, at half- filling, second- 
order perturbation theory in t, valid for small t/Uo, must 
yield an effective antiferromagnetic coupling proportional 
to t 2 /Uo- Thus, for small t/Uo, t 2 /Uo should equal the 
same constant at all antiferromagnetic phase transitions 
at half-filling (Recall that all of our coupling constants 
are dimensionless, incorporating the inverse temperature 
factor 1/kT). Equivalently, t/Uo should be linear in 1/t 
at all antiferromagnetic phase transitions at half-filling, 
for small t/Uo and therefore for small 1/t (low tempera- 
ture): 

1/t ~ t/Uo ■ (16) 

This is indeed rendered by our calculation, as seen in 
Fig. 3. For higher values of 1/t, Ea. (|16fl is not appli- 
cable, since second-order perturbation theory does not 
hold, and indeed our calculated curve in Fig. 3 deviates 
from linearity. On the other hand, the approximation 
in our recursion relation is even more justified, since the 
commutation relations that are ignored involve terms of 
order t 2 . 

The antiferromagnetic transition temperature as a 
function of coupling Uo/t is also shown in Fig. 4, together 
with calculated values from other approximation schemes 
for the d = 3 Hubbard model. We see that our results for 
intermediate coupling are comparable to those of quan- 
tum Monte Carlo studies |33Ll34| . As expected, our tran- 
sition temperature vanishes in the limit Uo/t — > 0, since 
there are no phase transitions for the non-interacting sys- 
tem. Thus, our approximation behaves correctly both at 
strong coupling (previous paragraph) and at weak cou- 
pling. 

THb and Ttj phases: For large values of Uo/t, the 
novel phase found in the tJ model 13] (which we call 
Ttj) also occurs in the Hubbard model. In addition, we 
find a closely related phase (Tjjb), unique to the Hubbard 
model, at smaller Uo/t. The two phases are characterized 
by very similar properties: the hopping strengths to, t 2 , 
and t x renormalize to ±oo, and the phase sinks have a 
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Chemical potential fj. /U 



FIG. 1: d — 3 Hubbard model phase diagrams in temperature versus chemical potential. The hole-rich disordered (hD), 
near- half-filled disordered (nHD), half-filled disordered (HD), electron-rich disordered (eD), antiferromagnetic (AF), Tub, and 
T t j phases are seen. The full curves are second-order phase boundaries, while the dotted curves are first-order boundaries. The 
dashed curves are not phase transitions, but disorder lines between the near-half-filled disordered and the hole-rich or electron- 
rich disordered phases. The progression (a) Uo/t = 20 through (e) Uo/t — 5 shows the changing phase diagram topology from 
strong to intermediate coupling. The r t j phase, which is prominent at strong coupling, disappears entirely for Uo/t < 6, and 
the THb phase is prominent for intermediate couplings. 
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FIG. 2: d = 3 Hubbard model phase diagrams in temperature versus electron density. The full curves are second-order phase 
boundaries. The coexistence boundaries of first-order transitions are drawn with dotted curves, with the unmarked areas 
inside corresponding to coexistence regions of the two phases at either side. The dashed curves are not phase transitions, but 
disorder lines between the near-half-filled disordered and the hole-rich or electron-rich disordered phases. Noteworthy is the 
narrowness of the first-order transitions, with jumps in the electron density of the order of a few percent (i.e., the width of the 
coexistence region). The antiferromagnetic phase is unstable to about 8-15% hole (or electron) doping away from half-filling. 
In the intermediate Uo/t regime, the r Hb phase appears for about 10-18% hole (or electron) doping. At larger Uo/t, the r tJ 
phase dominates, and exists between 30-35% hole (or electron) doping. 
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0.00 



0.05 



0.10 

t/U 



0.15 



0.20 



FIG. 3: The data points are the calculated antiferromag- 
netic transition temperatures at half-filling. The linear rela- 
tion that is expected for strong coupling at low temperatures 
(Sec.V) is obtained. 
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FIG. 4: Comparison of the antiferromagnetic transition tem- 
peratures at half-filling for the d — 3 Hubbard model cal- 
culated from various approaches: the renormalization-group 
method of the present paper (solid line); QMC [33| (dia- 
monds); QMC (triangles); DMFT [a| (squares); RPA 
weak-coupling expansion 36] (dot-dashed line); and two 
approximations for the strong-coupling behavior — the high- 
temperature expansion of the Heisenberg model, 1/t — 
3.83t/[/ HI (dotted line) and Weiss mean- field theory, 1/t — 
6t/Uo (dashed line). 



non-zero vacancy hopping expectation value 



-2/3 (r tJ ) 
0.663972 (r Hb ) 



(17) 

for Hq/Uq < 1/2, and a non-zero pair hopping expecta- 
tion value 



ni-a-Tlj-aic^Cja + c] CT C lCT )^ = i 



2/3 (r tJ ) 
-0.663972 (r Hb ) 



(18) 

for Ho/Uo > 1/2. In both cases, as expected for the 
occurrence of hopping, the electron densities at the sinks 
have values different from (empty), 1 (half filled), or 
2 (doubly occupied): (rij) = 2/3, 4/3 and (riiiij) = 1/3, 
5/3 for the r tJ phase, and (n t ) = 0.668014, 1.331986 and 



(mrij) = 0.336028, 1.663972 for the r Hb phase. (At the 
sinks of non-r phases, the electron density is, on the other 
hand, 0, 1, or 2.) There are also small spin correlations 
in the phase sink limits, (Si ■ Sj) = —1/4 for r t j, and 
(Si ■ Sj) = 0.0840069 for r Hb , which yield, throughout 
these phases, small antiferromagnetic correlations in the 
original system. 

The boundaries in Fig. 1 are controlled by fourteen 
unstable fixed points, given in Table VI. For smaller val- 
ues of Uo/t, the topology of the phase diagram is that 
of Fig. 1(e), where the AF/HD, AF/nHD, AF/r Hb , and 
hD/rnb boundaries are respectively controlled by the 
second-order fixed points C*, C^, C?j, and CJ. The lat- 
ter three boundaries intersect at the multicritical point 
B2, controlled by the fixed point Bj. A segment of the 
hD/nHD boundary just above this intersection is second- 
order, controlled by the fixed point C5, ending at the 
multicritical point Bi, controlled by the fixed point B*. 
The high-temperature section of the hD / nHD boundary 
is a disorder line, controlled by the null fixed point N*, 
i.e., there is no phase transition above Bi. 

As Uo/t is increased, the phase diagram topology be- 
comes more complex. For Uo/t > 6, the Ttj phase ap- 
pears, its boundary with hD controlled by the second- 
order fixed point Cg . Portions of the lower-temperature 
boundary between the hD and nHD phases become first- 
order (fixed point F*), and islands of AF appear above 
the Ttj phase; their boundaries with hD are also first- 
order (fixed point FJj). The intersections of these first- 
order boundaries with other phase boundaries are con- 
trolled by the additional multicritical points and B\ , 
and by the critical endpoint L* |3l| . 

As the coupling Uo/t varies, a most interesting aspect 
of the changing phase diagram topology is the relative 
sizes of the ThB and r t j phases. The ThB phase is largest 
at intermediate values of Uo/t, and gradually decreases in 
size as we move into the strong-coupling regime, break- 
ing up into narrow slivers until at large values of Uo/t 
only tiny remnants of it are left in the phase diagram. 
The Ttj phase appears at intermediate values of Uo/t, 
grows in size as Uo/t is increased, and occupies a promi- 
nent place in the diagram next to the AF phase in the 
strong-coupling regime. As discussed in Section VII, this 
is precisely what we expect, since the Hubbard phase 
diagram should approximately reproduce the t J model 
results ^3| in the large Uo/t limit. 

Phase diagrams in terms of temperature versus elec- 
tron density (rij) are shown in Fig. 2. It is seen that the 
antiferromagnetic phase is unstable to at most 15% hole 
(or electron) doping at low temperatures. The r Hb and r tJ 
phases exist at different doping values, with r Hb appear- 
ing for approximately 10-18% doping, directly adjacent 
to the AF phase, and t u in the 30-35% doping range. 
The narrowness of the first-order transitions, with jumps 
in the electron density of the order of a few percent, is 
noteworthy. 



11 



Phase sink 


Interaction constants 


Additional 
properties 


to 


tl 


*2 








J 


V2 


V 3 


v 4 


hole-rich 
disordered hD 














oo 


— oo 

H/U = const. 
















near-half- filled 
disordered nHD 


2 In 3 





-2 In 3 


oo 

fa 0.24(7 


oo 


oo 

fa 0.62(7 





— oo 

fa -0.47(7 


OO 
S3 0.69(7 


-oo 

fa -1.38(7 


(7-2^-2 V" 2 — V 3 = 
2V 3 +V 4 =0 


half-filled 
disordered HD 














oo 


OO 

= if 














(7-2^=0 


electron-rich 
disordered eD 














oo 


oo 

jU/t/ =const. 
















ant i ferro- 
magnetic AF 

Cmo/^o # 1 / 2 ) 


— oo 

S3 -0.29J7 





OO 
S3 Q.29C7 


oo 

fa 0.14(7 


oo 


OO 
M 0.57(7 


oo 

fa 0.29(7 


— oo 

fa -0.071(7 


OO 

v 3 /v — o 


— OO 
v 4 /u -» 


[7— 2^i— 2V 2 — V 3 — >0 
2V3+V4-.O 


ant i ferro- 
magnetic AF 

(vo/Uq = 1/2) 














oo 


OO 

=3ia 


OO 

J/U — * 


— OO 


oo 

w Jj 


— OO 


[7— 2//-2V 2 — V 3 =0 
2V 3 +V 4 =0 
*2-*0=0 


Hb 

(V0/U < 1/2) 


— oo 





OO 


— oo 


oo 


oo 

-if 


oo 

S3 It, 


— OO 

--if 


— oo 

V 3 /(7 — 


oo 

V4/E7 


as-4.35 


Hb 

(n/U > 1/2) 


— oo 





oo 

-iff 


— OO 

--if 


oo 


oo 

fa (7 


oo 

-if 


— oo 

--if 


OO 

V 3 /(7 — 


OO 

v 4 /u -» 


-t 2 +U-n+^J 
-V 2 —V^—^.35 


T t.I 

(CO/UO < 1/2) 


oo 

S3 0.131/ 





— oo 

fa -1.46(7 


oo 

0.52(7 


oo 


— OO 
f» -0.022(7 


— oo 

f» -0.87(7 


— oo 

S3 -0.50(7 


— OO 
fa -1.13(7 


— OO 
fa -0.21(7 




T t,I 

(fO/^0 > 1/2) 


OO 
S3 0.42U 





— oo 

ESS -0.038(7 


OO 
pa 0.15(7 


oo 


OO 
fSS 1.62(7 


— oo 

fa -0.25(7 


— OO 
S3 -0.861/ 


oo 

:a 0.39(7 


— OO 
:a -0.060(7 





Phase sink 










Expectation 


values 










(To) 


(Ti) 


(T2) 


(T*) 






(S z -S 3 ) 


(riiUj) 


(»>T™«1™j) 


(n«T n U n 3T n il) 


hD 
































nHD 

















1 





1 








HD 

















1 





1 








cD 














1 


2 





4 


2 


1 


AF 

















1 


1 

4 


1 








T Hb 
T tJ 


( 0.663 
I 972 






/ 

<, -0.66 
I 3972 

{i 






< 0.331 
I 986 

{1 


r 0.668 
1 014 
| 1.331 
I 986 

11 


0.0840 
069 

1 

4 


C 0.336 
1 028 
I 1.663 
I 972 

t 


4 0.331 
I 986 

{1 







TABLE V: Interaction constants and expectation values at the phase sink fixed points. For r Hb and r tJ , the values for 
Mo/(7o>| are given. The hopping expectation values (T a ) are: (To) = ^ CT (foi_ CT hj_ CT (4 CT c J(T +cj CT c 4(T )), (Ti) = I] CT ((n 4 _ CT 7ij _ CT + 
/i i _ CT n J _ CT )(c l t <J c JCT +ct .c« r )} ) (T a ) = £ CT (m_ CT n,; -"( c L c i* + C J CT C ^)>. < T *} = < C I T c j T c Ix c J-t + c] T c iT cj iCii ). 



VI. SPECIFIC HEAT RESULTS 

From the calculated expectation values of the oper- 
ators occurring in the Hamiltonian [Eq.pJ], we have 
obtained the dimensionlcss internal energy per bond 
(/3H(i,j)). Recall that dimensionless coupling constants 
are exhibited in the Hubbard Hamiltonian of Eq.Q, e.g., 



t 



(19) 



where t is a constant that does not depend on tempera- 
ture. The specific heat is calculated with 



C = 



d(H(i,j)) 
dT 



( d 

= k B I -Q—[{c\ a c, ja + c\ a c ia ) 

+ 7 dFi( n ^ niL + n ^ n n) I • ( 2 °) 



The partial derivatives are taken at fixed Uo/t and at 
fixed density (rij). 



In Fig. 5 we plot 7 = C/T for Uo/t = 15 at several 
different electron densities. (The corresponding phase 
diagram is shown in Fig. 2(b)). At half-filling, (m) 
1 .00. we observe a broad peak near the HD/AF transition 
temperature, which we can attribute to the onset of spin 
order. As we dope the system with holes, this peak gets 
sharper, becoming most pronounced near (rii) = 0.68, 
directly above the transition temperature between the hD 
and r tJ phases. In fact, the C/T curve shows a multipeak 
structure near the transition, a general characteristic of 
the phase diagram region just above the r tJ phase. At 
electron density {n/) = 0.60, no longer in the r tJ range, 
the peak decreases in size and broadens out again. 

The distinct nature of the r tJ and r Hb phases becomes 
clear when we look at the low temperature specific heat. 
In Fig. 6 we plot the coefficient 7 = C/T as a function 
of electron density for Uo/t = 15 and at low temperature 
1/t = 0.085. In the limit as T — > 0, 7 is a measure of 
the linear contribution to the specific heat due to quasi- 
particle excitations. Near half-filling, 7 is close to zero, 
increases to a small level with sufficient hole doping, falls 
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TABLE VI: Unstable fixed points. The fixed points of the fj,o/Uo < 1/2 half space are given here. 



to near zero again in the r Hb phase, and dramatically in- 
creases only after the system makes a narrow first-order 
transition to the hole-rich disordered phase. The steady 
rise of 7 in the hD phase with further hole doping is con- 
sistent with a Fermi liquid interpretation of this phase. 
The increase in 7 is interrupted by the r t J interval, where 
the curve makes a sharp oscillation, but continues in the 
hD region on the other side. 

We see that the r tJ phase has non-zero 7 at low tem- 
peratures, while the r Hb phase does not. In Figs. 7(a) 
and (b) we contrast the two r phases directly, comparing 
representative C /T curves for r tI and r Hb transitions. We 
observe that in the r Hb phase the low-temperature spe- 
cific heat exhibits an exponential form characteristic of 
a gap in the quasiparticle spectrum. Specific heat data 
points for temperatures l/t < 0.2, shown in the top right 
inset of Fig. 7(b), were found to fit a theoretical curve of 
the same form as in the T — > limit of a weakly-coupled, 
BCS-type superconductor, 

C A ( A\ 

fcj; = ^372" ex P \f) ' (21) 

with a best-fit coefficient A = 1.02 ± 0.06 and a zero- 
temperature gap A = 1.01±0.01, where t~ x is used as the 



temperature variable. In contrast, the r tJ phase clearly 
has a gapless spectrum, as we see in the C/T curve of 
Fig. 7(a). As mentioned earlier, we also clearly see mul- 
tiple peaks in the specific heat just above the hD/r tJ 
transition temperature. 

The r tJ and r Hb phases have similar properties at 
the phase sink, most notably a non-zero hopping am- 
plitude, and thus are both good candidates for super- 
conductivity. Since the two phases are dominant in dif- 
ferent Uo/t regimes, their contrasting specific heat char- 
acteristics can potentially be understood as the differ- 
ence between strongly-coupled and weakly-coupled su- 
perconducting phases. For the strongly- coupled, BEC- 
like case, pairing occurs above T c , and these tightly 
bound bosonic pairs condense at the transition temper- 
ature. The double-peak structure in the specific heat 
above the r tJ phase is a possible indicator of such pair 
formation. Additionally, we expect that a BEC-like su- 
perconducting transition in three dimensions should have 
a specific heat critical exponent a = — 1 24] . Analysis 
of the Cg fixed point, governing the hD/r tJ boundary 
yields the result a = —0.97. The presence of low-lying 
excitations in a Bose gas is also consistent with the fact 
that we do not see a gap in the low-temperature specific 
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FIG. 5: The specific heat coefficient 7 = C/T as a function of 
temperature for Uo/t = 15, at several different electron densi- 
ties (rii) indicated in the legend. For this temperature range 
the densities (rn) — 1.00 and 0.93 lie inside the antiferromag- 
netic (AF) phase, 0.84 inside the r Hb phase, 0.73 and 0.60 
inside the hole-rich disordered (hD) phase, and 0.68 inside 
the T tJ phase. Here and in the following figures, 7 is shown 
in units of k% /i, where i is the temperature-independent con- 
stant in Eq. fB^ . 
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Electron density <«,) 
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FIG. 6: The specific heat coefficient 7 = C/T for Uo/t = 15 
at the low temperature of l/t = 0.085, as a function of elec- 
tron density (rii). The corresponding phases are indicated 
near the top of the figure, with second-order phase bound- 
aries marked by thin vertical lines. The interval between the 
vertical dashed lines corresponds to the first-order phase tran- 
sition. 




0.2 0.3 
Temperature 



FIG. 7: The specific heat coefficient 7 = C/T as a function 
of temperature for two different electron densities and values 
of Uo/t: (a) (m) = 0.68, Uo/t = 20; (b) (n t ) = 0.875, Uo/t = 
7.5. Phases are indicated near the top of the figures, with 
second-order phase boundaries marked by thin vertical lines. 
The interval between the vertical dashed lines corresponds to 
a first-order phase transition. In diagram (b) the top left in- 
set shows a close-up of the cusp in 7 at the AF/r Hb transition 
temperature. The data points in the top right inset are calcu- 
lated 7 values for temperatures l/t < 0.2, fitted to a BCS-like 
exponential curve of the form C/fesT = ^ 2 exp ( — -fr) , with 
best-fit parameters A = 1.02 ± 0.06 and" A = 1.01 ± 0.01, 
where t' 1 is used as the temperature variable. 



[37| | (examples of transitions in this class include the su- 
pcrfmid transition of 4 He, the superconducting transi- 
tion in certain high-T c materials like Y-123, and also in 
conventional superconductors, though for the latter the 
critical region is too narrow to be observed experimen- 
tally) [24|. Our calculated a is closer to the d = 3 XY 
than to the BEC value, supporting the weak-coupling 
interpretation of the r Hb phase. 



heat of the r tJ phase. 

Turning now to the r Hb phase, we already noted that 
its specific heat can be closely fitted at low tempera- 
tures to a BCS-like exponential curve, which is exactly 
what we would expect for a weakly-coupled supercon- 
ducting phase. Analysis of the C3 fixed point, control- 
ling the AF/r Hb boundary, yields a specific heat coeffi- 
cient a = —0.27. This translates into a finite cusp at the 
transition temperature, as shown in the top left inset of 
Fig. 7(b). For weak and intermediate couplings the su- 
perconducting transition is expected to belong to the uni- 
versality class of the d = 3 XY model, with a = —0.013 



VII. THE tJ LIMIT OF THE HUBBARD MODEL 

In the strong-coupling limit Uq 3> t, second-order per- 
turbation theory in t/Uo applied to the Hubbard model 
leads to the followin g H amiltonian (known as the tJ 

model) [n m m miii in m m hi, 

+ J^2(S i -S j -^n i n j ), (22) 

(ij) 
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FIG. 8: d = 3 Hubbard model phase diagram for large 
Coulomb repulsion Uo/t = 50 in temperature versus (a) chem- 
ical potential, (b) electron density (n;). The full curves are 
second-order phase boundaries, while the dotted curves indi- 
cate first-order boundaries. The dashed lines are not phase 
transitions, but disorder lines between the near-half-filled dis- 
ordered and hole-rich disordered phases. 




0.2 

Temperature l/t 

FIG. 9: The specific heat coefficient 7 = C/T as a function 
of temperature for Uo/t = 50 and (rij) = 0.67. Phases are 
indicated near the top of the figure, with second-order phase 
boundaries marked by thin vertical lines. The dashed lines 
are not phase transitions, but disorder lines between the near- 
half-filled disordered and hole-rich disordered phases. 
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Electron density (n.) 

10: The kinetic energy per bond (K) = — 



1.00 



+ 



} as a function of electron density (m) at temperature 
i/t = 0.2 for Coulomb repulsions U /t = 20, 50, 100, and 
1000 (indicated by numbers next to each curve). The solid 
curve at the bottom is the result calculated using the t J model 
renormalization-group equations |13lll4|| at the same temper- 
ature, with the corresponding J/t — 0.004. 



where J = At 2 /U and P is a projection operator prohibit- 
ing double occupation of a lattice site. In addition to the 
terms shown above, the perturbation theory generates 



a three-site term of the form 



(ikj) 



(Sfc)c 



but 



this term is usually ignored, from the assumption that 
it does not radically alter the physics of the tJ model. 
(Our current results, directly from the strong-coupling 
limit of the actual Hubbard model, confirm this assump- 
tion.) We thus expect that our Hubbard model approach 
in the limit of large Uo/t should give results qualita- 
tively similar to those found for the tJ model in earlier 
renormalization-group studies [HEI. The phases of the 
tJ model found in these studies are identical to those of 
the Hubbard model, except that there is no r Hb phase. 

Fig. 8 shows the Hubbard model phase diagram in 
terms of temperature versus chemical potential and tem- 
perature versus electron density for Uo/t — 50. At this 
large coupling, we do indeed observe a phase diagram 
very similar to that found in the earlier study of the t J 
model m>mi- I n particular, the r t j phase is surrounded 
by AF islands, and directly above Ttj we get a lamellar 
structure of alternating AF, nHD, and hD phases. The 
AF phase near half-filling is unstable to only about 5% 
hole doping. This phase diagram can be seen as an evo- 
lution from the Uo/t — 20 result of Figs. 1(a) and 2(a), 
with the ThB entirely disappearing at Uo/t = 50 except 
for infinitesimal slivers. The multiple peaks in the spe- 
cific heat above the r t j transition persist in the strong- 
coupling limit, as seen in Fig. 9, which plots the specific 
heat coefficient 7 for U /t = 50 at {n/) — 0.67. The peak 
structure here is more complex than in Fig. 7(a), due to 
the above-mentioned lamellar phases. 

We can also observe the evolution from the Hubbard to 
the t J limits through the expectation value of the kinetic 



energy per bond, (K) = - Y,A c ia c : 
proportional to the density of free carriers in the system. 



c] a c ia ), which is 
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Fig. 10 shows (K) as a function of electron density for 
the temperature 1/t = 0.2, calculated at several different 
couplings Uo/t. As Uo/t is increased, the value of (K) 
at half- filling is reduced, and when Uo/t = 1000 we are 
close to the t J limit, with the kinetic energy at half-filling 
almost zero, indicating no available free carriers due to 
the prohibitively high energy of double occupation. The 
Uo/t = 1000 curve almost exactly overlaps the result cal- 
culated from the t J model renormalization-group equa- 
tions at the same temperature using the corresponding 
coupling J/t = At/Uo = 0.004. 
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APPENDIX A: DETERMINATION OF THE 7p 
IN TERMS OF THE MATRIX ELEMENTS 
OF THE THREE-SITE HAMILTONIAN 

Eq.(8) allows us to express the matrix elements "f p = 
(</> p |e~ /3 H ( l,fc ) \(j> p ) of the renormalized, exponentiated 
two-site Hamiltonian in terms of matrix elements of the 



unrenormalized, exponentiated three-site Hamiltonian, 
as given below. The 7 P , in turn, determine the renormal- 
ized interaction constants, in Eq.(10). In the equations 
below, (^ q \\tp q ) denotes (^ q \e-^ H ^-^ H ^\^ q ): 

71 = (0l||0l) + 2(02||02) + (0 9 ||09), 

72 = (03 1 I V>3> + \ (08 1 108 ) + \ (015 I |01 5 ) + (024 1 |02 4 ), 

1 3 

74 =(0 6 | 106 ) + ^ (0n I I01l) + 2 (^20 |02o) + (032 life), 

76 = (0io||0io) + 2(0 26 | |02e) + ("0441 |044>, 

77 =<"013| 1013) + (034 1 1^34) + <V'38|iV'38> + (048 I |04s), 

78 = (012 1 1012) + 2(0 3i 1 |03l) + (^47 1 1047), 

2 4 

79 = (014 1 1^14) + jj ("023 I + g (039 I |03g) + (049 | |04g) , 

1 3 

7l 2 = (025 1 1 025 ) + 75 (045 1 1 045 ) + 77 (050 1 1 050 ) + (059 1 1 059 ) , 

1 3 

714 = (033 1 1 033 )+ 2 (046 1 1 046 } + - (055 1 1 055 ) + (062 1 1 062 } , 
716 = (043 1 1 043 ) + 2 (058 1 1 058 } + (064 1 1 064 ) , 

= (01O| 1012) +2(0 26 ||03l} + (044H047)- 
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